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ABSTRACT 


A 2-D finite element model is presented for the melt growth of single crystals in a 
microgravity environment with a superimposed DC magnetic field. The model is developed 
based on the deforming finite element methodology and is capable of predicting the phenomena 
of the steady and transient convective flows, heat transfer, solute distribution, and solid-liquid 
interface morphology associated with the melt growth of single crystals in microgravity with and 
without an applied magnetic field. Numerical simulations were carried out for a wide range of 
parameters including idealized microgravity conditions, the synthesized g-jitter and the real g- 
jitter data taken by on-board accelerometers during space flights. The results reveal that the time 
varying g-jitter disturbances, although small in magnitude, cause an appreciable convective flow 
in the liquid pool, which in turn produces detrimental effects during the space processing of 
single crystal growth. An applied magnetic field of appropriate strength, superimposed on 
microgravity, can be very effective in suppressing the deleterious effects resulting from the g- 
jitter disturbances. 


NOMEMCLATURE 


U : velocity vector 

u : x direction velocity component 

v : y direction velocity component 

p : pressure 

T : temperature 

C : concentration 

T m0 : melting temperature of pure Sn 

p 0 : density of pure Sn at T m 

C p : specific heat of Sn 

v : kinematic viscosity 
D : solute diffusivity 
j5 T : thermal expansion coefficient 
[3 C : solutal expansion coefficient 
k 0 : segregation coefficient 
m : slope of solidus curve 
g g : earth gravity constant 
H : latent heat of Sn 
k : thermal conductivity 
k : thermal diffusivity 
s : black body coefficient of ampoule 
W t : inner width of ampoule 
W 0 : outer width of ampoule 
L : length of three zones 
L h : length of hot zone 

L c : length of cold zone 

L g : length of gradient zone 

dT/dx : temperature gradient 
T h : hot zone temperature 

T c : cold zone temperature 

V g : crystal pulling speed 

B 0 : magnetic field magnitude 

U 0 : velocity scale 

A T : temperature scale 

R 0 : length scale 

t 0 : time scale 


C 0 : concentration scale 

h{y) : dimensionless x - position of 

growth interface 

L 0 : dimensionless x - position of hot 
end 

L, ,L t : dimensionless lower and upper y - 
positions of ampoule outer surface 
L 2 , 1 3 : dimensionless lower and upper y - 

position of ampoule inner surface 
Pr : Prandtl number 

Sc : Schmidt number 

Gr T : thermal Grashof number 

Gr s : solutalGrashof number 

Ste : Stefan number 

Ha : Hartmann number 

Ra : radiation number of ampoule 

T m rdimensionless melting point 

V p .'dimensionless pulling speed 

O' : interface velocity 

h :unit normal vector 

i,j : unit vector of i th, j th component 
s : penalty parameter for pressure 

Q : computational domain 

dQ : computational domain boundary 

Subscripts 

cry : crystal 

amp : ampoule 

amb : ambient 

i,j : /th,/'th component 

Superscipts 

T : matrix transpose 



1. INTRODUCTION 

Space flight experiments on the melt growth of single crystals indicate that natural 
convection induced by gravity perturbation or g-jitter has a strong influence on the defect 
formation in the crystals. One of the important characteristics of g-jitter-induced convection is 
that it occurs randomly in both direction and time, which is attributed to the nature of the g-jitter 
forces. During normal space flights, a gravity level of 10' 6 g Q is established. This dynamic 
weightlessness condition experienced by flight vehicles is prone to perturbations that stem from 
a variety of sources, which include the crew motions, mechanical vibrations (pumps, motors, 
excitations of natural frequencies of spacecraft structures), spacecraft maneuvers and attitude, 
atmospheric drag and the Earth's gravity gradient [1]. Detailed studies on g-jitter effects reveal 
that g-jitter induced convection in spacecraft is directly related to the magnitude and frequency 
of the gravity perturbation and to the alignment of the gravity field with respect to the growth 
direction or the direction of the temperature gradient [2-6]. Crystal growth experiments further 

indicate that a residual gravity of 1 O' 5 g 0 is sufficient to cause appreciable fluid motion in the 
liquid, thereby resulting in the quality of the crystals being unacceptable. The orientation of the 
gravity vector with respect to the temperature gradient also plays an important role in convective 
melt flows. The velocity attains a maximum when the gravity vector is perpendicular to the 
temperature gradients. 

Numerous numerical simulations have been conducted to understand and estimate the 
adverse effer.ts of time, varying g-jitter [6-10]. Both 2D and 3D numerical models have been 
developed for this purpose [4,5,1 1]. These models have been used to study the effects associated 
with both idealized single- and multiple-frequency g-jitter modulations and real g-jitter data 
collected by accelerometer during actual flight experiments. These studies showed that the 
frequency, amplitude and spatial orientation of the residual gravity vector all play an important 
role in determining the convective flow behavior of the system. When the residual accelerations 
oscillate about the positive and negative of an axis, the orientation of this direction relative to the 



density gradient determines whether a mean flow is generated in the system [12]. Sinusoidal 
accelerations induce an oscillating convective flow and composition oscillation in the liquid. 
The resulting velocity and composition fields oscillate with the same frequency of the affecting 
gravity field. Benjapiyapom et al. [13] reported a study on the g-jitter influence of solute 
transport and convection together with solidification in a planned space flight experiment 
configuration. They used the finite difference approximation for field calculations and the 
enthalpy method to model the solid-liquid interface. In their calculations, the melting 
temperature of the pure solvent was used to determine the solid-liquid interface. While the 
enthalpy method is very useful for alloy solidification with a distinctive mushy zone, it is in 
general considered inadequate for precisely capturing the solidification front of pure or nearly 
pure materials. They used simulated g-jitter data for their simulations to show that for large 
frequencies, the higher amplitude of the gravitational acceleration is required to produce effects 
on the segregation. Their study further illustrates that the concentration varies appreciable along 
the solid-liquid interface, suggesting that it is necessary to take into consideration the effects of 
concentration on the melting temperature to get an accurate estimate of the interface shape. 

The use of the magnetic field to control the melt motion during single crystal growth has 
been widespread in the semiconductor industry. Because the magnetic damping effect comes 
from a different origin from that of gravity reduction, it is plausible to suppress the g-jitter 
effects by imposing an appropriate magnetic field [14]. Studies on the g-jitter-induced melt 
flows in the systems of relevance to space crystal growth experiments have been reported by Ma 
and Walker [12] and Baumgartl and Muller [11, 15]. The former presented a semi-analytical 
model for the magnetic field effects on g-jitter driven thermal convection, while the latter 
presented a 3-D model of melt flow driven by a single frequency g-jitter force parallel to the 
thermal gradient with either liansveise or axial fields imposed. Recently, Li and his co-workers 
[7,14,16-20] have reported a series of studies, including both analytical analyses and numerical 
simulations, on the magnetic field effects on the g-jitter induced flows and mass transfer. Unlike 
others, their studies used both simulated g-jitter signatures and real g-jitter data taken from 
typical space flights, with the latter intended to obtain an understanding of the thermal and flow 
phenomena in a realistic g-jitter environment. In space vehicles, g-jitter is not aligned perfectly 
with the thermal gradient during crystal growth. The most deleterious flow effects, however, 


result from the g-jitter component that is perpendicular to the thermal gradient [4]. The melt 
flow in microgravity is shown to have a very strong effect on the solute transport. Their studies 
further demonstrate that these deleterious g-jitter effects can be suppressed by the application of 
an external magnetic field. Thus far, these studies have been primarily concerned about the 
convection with or without solute striation in the presence of magnetic fields in the melt pool for 
a configuration, which is a version idealized, though relevant, version of a melt growth 
experiment furnace used for space flights. None of these studies have considered the 
solidification behavior, which is a critical part of the process. Nonetheless, these studies have 
provided a useful basis upon which models incorporating more realistic issues such as 
solidification and ampoule configuration, coupled with more realistic thermal boundary 
conditions. These more realistic models will provide a much more powerful tool not only for 
enhancing our fundamental understanding of fluid mechanics and thermal phenomena in space 
crystal growth systems but also for developing rational guidelines for both space experiment 
design and interpretation of experimental measurements. 

In the present study, a finite element model is presented for the heat and mass transfer and 
solidification phenomena associated with the melt growth of single crystals in microgravity with 
and without an imposed magnetic field. The analysis is based on a 2-D geometry. It is noted 
that a furnace for the melt growth is of cylindrical geometry and under terrestrial conditions a 2- 
D axisymmetric condition can be established and maintained for field variables of interest [21]. 
In a microgravity environment, however, such an axisymmetry condition can not be maintained 
because the frequent maneuvering of space vehicles changes the gravity orientation constantly. It 
is also realized that because of the frequent change in gravity directions, in particular during the 
period when g-jitter sets in, the flow structure is truly three-dimensional and 3-D models should 
be developed to fully appreciate the complex 3-D flow structure [21]. Nonetheless, a 2-D 
analysis will provide some essential features of g-jitter induced flows and solidification behavior 
and thus form the required basis upon which 3-D models can be further developed. The 2-D 
analysis presented below is based on the solution of the Navier-Stokes equations for fluid flow, 
the energy balance equation for the temperature distribution and the mass balance equation for 
species transport in an Sn-doped Bi alloy melt. The numerical scheme entails the 
implementation of a quasi-Lagrangian deforming finite element approach where the solid-liquid 


interface is tracked precisely by deforming the elements associated with the moving interface. 
This interface track approach allows us to accurately predict the interface shape change as a 
function of solute distribution, which is critical for space experiments but has not yet fully 
appreciated. Numerical simulations are conducted of the transient convective flows and solute 
redistribution in the melts and their effects on the solidification interface morphology under the 
combined action of g-jitter and external magnetic fields. Both synthetic g-jitter, which 
represents a single component of Fourier synthesis of g-jitter data, and the real g-jitter data are 
used in the analyses. 

2. PROBLEM STATEMENT 

Figure 1 schematically illustrates the 2-D model used to study the Bridgman-Stockbarger Sn- 
doped Bi single crystal growth in microgravity. An external DC magnetic field is applied along 
the thermal gradient. The ambient temperature profile ( T amb ) consists a cold zone (T e ), a linear 
gradient zone ( dT/dx ) and a hot zone (T h ). In the present study, the pseudo-steady-state model 
(PSSM) [22] is adopted in which the translation of the ampoule is simulated by feeding the melt 
with the uniform impurity concentration C 0 into the inlet at the constant growth velocity V g , 
while the crystal is withdrawn from the outlet at the speed that conserves the mass of the system. 
The length of the computation domain is L , the inner width of the ampoule is W. t and the outer 

width is W 0 . The configuration described above is similar to that used for the MEPHISTO-2 and 
MEPHISTO— 4 Space Flight experiments [6,23]. 

The residual acceleration experienced by the spacecraft in microgravity is composed of a 
steady acceleration component and a time-dependent acceleration component (g-jitter). The 
amplitude of the time-dependent component is much greater than that of the steady component. 
The real g-jitter data taken during space flight show that g-jitter changes both spatially and 
temporally. In some cases, to help analyses, g-jitter can be represented by the combination of a 
series of single frequency accelerations as follow: 

N 

g(t)= Tin sin(2 7ico n t) 

77= 1 


( 1 ) 



where g n is the amplitude and co n the frequency associated with the nth component of g-jitter. 
Both the synthetic and real g-jitter data are considered in the present study. 

The time dependent phenomena of fluid flow, heat and mass transfer in the above system 
with the presence of an external magnetic field are described by the Navier-Stokes equations 
with the time varying gravity force and the Lorentz force resulting from the imposition of a 
magnetic field. For the melt flow, the standard Boussinesq approximation, 
p = p Q [l - P t (T - T m0 ) - /? c (C - C 0 )] , has been used. Furthermore, the thermoelectric effect, 

i.e., J = -SVT , is neglected, because of the lack of reliable data for the absolute thermoelectric 
power S of the fluid and also because it generates a flow component in the z-direction and thus 
violates the 2-D assumption. With an appropriate choice of scales, the governing equations may 
be non-dimensionalized and are given below. For the melt, these equations become 


v-u = o 

— + U-VU = -Vp + V 2 U-[Gr r (T-T m ) + Gr s (C-\)\g(t) + Ha 2 (UxBxB) 


- + U-VT = —V 2 T 

dt Pr 

- + U-VC = —V 2 C 

dt Sc 


( 2 ) 

(3) 

(4) 

(5) 


For the crystal and the ampoule, where a bulk flow is superimposed because of the choice of the 
coordinate system, only heat transfer analysis is needed. The solute diffusion in crystal is also 
ignored. With these simplifications, the governing equations can be written as follows, 


a Tcry + v P (e x • vr ) = V 2 Z 


d t 


c,y ' /r Pr 


cry 


dT 


amp 


+ Vp(e x -VT) = 


“amp 1 


V 2 T 


( 6 ) 


dt 


k Pr 


(7) 



In deriving the above equations, the following scales are used: R 0 = 0.5 IF, for length, 
U 0 = v/R 0 for velocity, t 0 = R 0 2 /v for time, p 0 = p 0 U 0 2 for pressure, C 0 for concentration, 
g 0 for gravity, B 0 for magnetic field, AT = (T h - T c ) for temperature while the temperature is 
non-dimensioned by (T - T C )/AT . The non-dimension parameters are defined as: 

Gr T = f3 T R 0 g 0 AT /v , Gv s = fl s R 0 g 0 C 0 /v , 

Pr = v/*r, Sc = vjD , Xa = s„, p aR t (ATf /k^, 

V, = VJV \ , Ste = H/ATC , , Ha = B^TTIpT . 

The boundary conditions for the configuration considered are determined by the physical 
constraints and are detailed below. 


T = 0, 

at x = 0, Zj < y < Z 4 

(8) 

T = 1, 

at x = Z 0 , Z, < y < Z 4 

(9) 


at 0 < x < Z 0 y = Z, and y = L A 

(10) 

n-VT = (k m , p /kJ(n-VT„ p ), 

at 0 < x < ft(y), y = Z 2 arccf y = Z 3 

(11) 

u = V p , v = 0, n ■ VC = 0 , 

along the melt-ampoule interface 

(12) 

«-vr = (w*)(«.vr.J. 

at /z(y) < x < Z 0 , y = Z 2 ancf y = Z 3 

(13) 

u — V p , v = 0 , h • VC = V p Sc(C - \)(n • 

), at x = Z 0 , Z 2 < y < Z 3 

(14) 

), « - (c 7 - V, 

)"0, at x - h(y), Z 2 < y < Z 3 

(15) 

n-VC = Sc[n-(i'-C?')c-«-(l7 w -C?' 

JC co ,], at x = /z(y), Z 2 < y < Z 3 

(16) 

T = T m +mC , 

at x = /z(y), Z 2 < y < Z 3 

(17) 

(a ■ vr) - (» • vr t „ ) = s«p r « • (t?,. 

,-£/*), at x = /z(y), Z 2 <y<Z 3 

(18) 


where T* = T C /A7\ C cry =k 0 C , n-U cry = V p> hxU cry =0 . 



3. THE FINITE ELEMENT SOLUTION 

The governing equations described above along with the boundary conditions are solved 
using the deforming Galerkin finite element method. The stiffness matrix is obtained by using 
Galerkin’s method of Weighted Residuals. The formulations and relevant benchmark tests were 
detailed in a series of papers published earlier [17, 24] and thus only a brief summary is given 
here. The governing equations are recast in an integral form and the field variables are 
interpolated using shape functions over the computational domain. With an appropriate algebraic 
manipulation, the following set of equations are obtained, 

(£y/f • V<f> T dV)V ,-=-*(£ W T dV)Y (19) 

/7f T 

( £ (fxj) T dV ) — *- + ( £ <f> u • Shj> T dV)V i - ( £ / • S/(f>Y T dV) P + ( £ • V<p T dV)Vj 

+ (£(?• • V^ r W)Vj + ( l<t>e T gdV)(T - T m ) + ( £ f w T gdV){c - 1) 

= \ x n-r-i<f>dS+ (20) 

(^Pr00 r dV)^- + (^PT0u-V0 T dV)T + (^V0-V0 r dV)T=- ^q^dS (21) 

( | Sc00 T dV) ^- + (^Sc0 u V0 T dV)C + ( £ V0- V0 T dV)C = - £ q c 0ds (22) 

Once the form of shape functions <j>, 0, and \|/ for velocity, pressure and scalars are specified, 
the integrals defined in the above equations can be expressed in matrix form. Combining the 
momentum and energy equations into a single matrix equation gives rise to the following 
element stiffness matrix equation, 


'M 

0 

0 ‘ 

if 


’aOUHK + JtEM^E 7- 

B r 

B c 

U‘ 

0 

N r 

0 

T 

+ 

0 

Dy (U) + L j 

0 

T 

_0 

0 

N c _ 

C 


0 

0 

D C (U) + L C 

C 


F 

G r 

G c 


(23) 


Note that in constructing the above element matrix equation, the penalty formulation has been 
applied, and P in the momentum equation is substituted by^M" 1 E r U. The assembled global 



matrix equations are stored in the skyline form and solved using the Gaussian elimination 
method. The coefficient matrices of Eq. (23) above are calculated by 


M, = j '^W T dV ; N r = JPr 6d T dV 

M =^ed T dV; N c = ^Sc00 T dV 

E, = • V<f>y/ T dV\ L c = ]V0 ■ Vd T dV 

L T = fV0-V0 T dV; A(U) = u • V 0 T dV 

D C (U)= [ i Sceu-V0 T dV; D r (U) = _[Pr0u- V0 T dV 

Bc^lGr^g^W; B r = £ Gr T (g </>9 T )dV 

G c =-^q c BdS\ G T =-^q T 6dS 

F = ^n-r<f)dS + [ i <f>Ha 2 {yL^xB)dV + ^GrcigtfiB^T^V + [ i Gr T (g<f>9 T )dV 
K, =([v < j>-VfdV)5 il + l (i-V0)(j-V0 T )dV 

To implement the deforming finite elements to model the dynamic change of the moving 
interface, i.e., solidification front between the liquid and solid, a quasi-Lagrangian description is 
adopted. By this method, a region that covers the solidifying liquid and solid is defined and the 
nodes within the region are allowed to move in accordance with the interface movement. These 
additional velocities that result from the mesh movement are added to the velocity field as given 
in the above equations. The energy balance equation describing the latent heat release and 
interface change (Eq. 1 8) is directly integrated within the context of weighted residuals, 

(J, en-V0 T dV)T-(l(k^/k)eh-Vff T dV)J ^ =- BdS (24) 

which is applied as a surface energy source to the total thermal energy balance equation, and 
added to the right-hand side of Eq. (23). Two basic algorithms have been applied in the past. 
One treats the moving interface as a separate variable and the interface thermal boundary 
condition in this present case would then be incorporated into the global finite element matrix. 



The other treats the interface thermal boundary condition as a separate constraint and solved 
separately from the coupled field equations. Our numerical experience with these types of 
moving interface problems and solution algorithms shows that in a majority of cases the 
incorporation of the isothermal constraint into the global matrix results in a very unstable 
nonlinear system. The whole system is extremely sensitive to the boundary perturbations, more 
often than not, leading to divergence. Separation of the moving interface boundary coordinates 
from the global finite element solutions for field variables, however, requires the assurance of 
convergence of both moving interface coordinates and field variables in two related loops. For 
this case the global finite element matrix has a smaller bandwidth and thus requires a short CPU 
time for each iteration of field variable computations, which takes the majority of the computing 
time. Our tests further showed that the two methods require approximately the same CPU time, 
with the latter being much more robust and numerically stable. The latter has been used in the 
present study and for the results presented below. Also, there are two different methods used to 
track the liquid-solid interface and the results are the same within the machine accuracy. One 
method entails the use of the weighted residual method to integrate the interface thermal 
boundary condition (Eq. (17)), while the other the use of a directional search of the melting 
temperature after the coupled thermal, solute and fluid flow field calculations are converged. 
The latter is adopted in the present study. 

In summary, our interface tracking strategy used in the present study involves an iterative 
procedure that entails applying the energy balance equation along the interface as a surface 
source and searching for the interface position coordinates based on each converged field 
calculations. The updated interface positions are the fed back to the field calculations until both 
the interface position coordinates and field variables are converged within a preset criterion, 
which is set at 1 x 10' 4 fur the results presented below. 

4. RESULTS AND DISCUSSION 

The finite element model described above permits the prediction of both steady and transient 
fluid flow, heat and mass transfer, along with the solid-liquid interface morphology in single 
crystal growth systems under both terrestrial and microgravity conditions with and without the 
presence of magnetic field in the axial direction. The model development was based on the finite 



element code that has been reported in early studies, for example [7,14,16-20], which include 
both code development and benchmark testing. While both the formulations and computer code 
are general and can be applied to both 2-D and 3-D calculations, our focus here is on a 2-D 
simplification of a cylindrical system for growing binary single crystals in a microgravity 
environment. The thermo-physical properties and geometric dimensions for the study are given 
in Table 1 . The studies used 2,968 9-node elements with the increasing mesh distribution near 
the interface. The penalty pressure formulation was used to approximate the pressure field. Tests 
were conducted to ensure the numerical results are mesh independent to within a tolerance of 
less 1% between two consecutively refined meshes. Numerical simulations were carried out for 
a wide range of conditions and a selection of the computed results is presented below. 

1> Steady microgravity condition 

Under an idealized microgravity condition, the gravity level is 10“ 6 g 0 . This condition may 
be obtained when the space vehicle follows the orbit perfectly without disturbances. Numerical 
simulations were carried out for this condition without an imposed magnetic field. The 
calculations will serve a dual purpose. The results are useful in helping to gain physical insight 
into the steady state behavior of the system. They also provide the initial condition for the 
dynamic behavior of the system when g-jitter effects set in. The simulations were based on the 
direct numerical solution to the steady state governing equations by dropping out the time 
dependent terms in the model equations described in Section 2 while keeping all the mutually 
coupled terms. For this case, the gravity is assumed to be oriented perpendicular to the growth 
direction, which represents the worst scenario case. The computed results show that the flow 
field is dominated overwhelmingly by the growth velocity, and the thermal and solutal gradient 
induced convective flow recirculation is small. These results are consistent with those obtained 
earlier from a simplified melt flow system where solidification was not considered [14], 
confirming that these simplified models indeed provide a good approximation to the melt flows 
under more realistic growth conditions, as they are intended to. Temperature distributions in the 
growth system are largely controlled by thermal conduction, which is a direct result of small 
Prandtl number for the system being studied. The solute distribution, one other hand, appears to 
experience a slight distortion, which results from the small convection induced by the 
microgravity force. 



Detailed analysis shows that for this case, the solute distribution has a striation of about 
7%, measured by (C max - C min )/C merage along the solid-liquid interface. This concentration non- 
uniformity, according to Eq.(17), will affect the melting temperature and thus the system 
performance. Figure 2a compares the computed results of concentration and the normalized 
melting temperature distributions along the solid-liquid interface, and the solidification interface 
position for the system with and without taking into account the concentration effect on the 
melting temperature. Their differences are also showed in Figure 2b. It transpires that with the 
concentration effects considered, the melting temperature is no longer a constant, which results 
in a noticeable change in the solidification interface position. It is thus important to include the 
concentration effect on the melting temperature. 

2> Single frequency g-jitter 

Single frequency g-jitter represents an idealized perturbation to a mcirogravity 
environment. The Fourier analysis of g-jitter data collected by on-board accelerometers indicates 
that g-jitter has a wide range of frequency spectrum. However, space-flight experiments with 
single crystal growth suggest that only the low frequency g-jitter produces deleterious effects on 
the quality of the crystals. A single component g-jitter force provides a valuable means for 
developing a fundamental understanding of physics governing the effects of gravity perturbation 
on convection and solidification phenomena associated with space melt growth of crystals. The 
computed results for the dynamic development of fluid field, temperature and solute distribution 
and solidification interface morphology are given in Figures 3 to 4 for a quasi-steady state time 

harmonic cycle. The results were obtained a single frequency g-jitter, g(t ) = 10~ 3 g- 0 sin(0.2/rt) , 
oriented in the direction perpendicular to the thermal gradient without an imposed magnetic 
field. For this case, a strong convection recirculation develops in the melt pool, which covers the 
entire pool region. The convection loop oscillates and reverses its rotating direction with time at 
approximately the same frequency as the applied g-jitter. A very strong flow sweeps along the 
solid-liquid interface. Note that the ampoule and the solid move at the crystal pulling velocity 
with respect to the chosen frame of reference, in contrast to the recirculating convective flow 
motion in the liquid (see Figure 3). The effect of the convection on the concentration 
distribution in the melt is shown in Figure 4, where an appreciable concentration striation is 



5. CONCLUDING REMARKS 

This paper has presented a 2-D finite element model for the melt growth of single crystals in 
a microgravity environment with a superimposed magnetic field. The model is developed based 
on the deforming finite element methodology and is capable of predicting the phenomena of 
steady state and transient melt convective flows, heat transfer and solute distribution, and solid- 
liquid interface morphology associated with single crystal growth in microgravity with and 
without an applied magnetic field. The field variables are solved in a completely coupled 
fashion, with allowance made for the effects of solute striation on the melting temperature. 
Numerical simulations were carried out for a wide range of parameters including the idealized 
microgravity condition, the synthesized g-jitter and also real g-jitter data taken by on-board 
accelerometers during a space flight. The results show that the time varying g-jitter disturbances, 
although small in magnitude, indeed cause an appreciable convective flow in the liquid pool, 
which in turn produces detrimental effects during space growth of single crystals. An applied 
magnetic field of appropriate strength, superimposed on microgravity, can provide a viable 
means to suppress or eliminate (i.e. suppress beyond detection) the deleterious effects resulting 
from the g-jitter disturbance. 
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Table 1. Parameters used for calculations 


Physics properties 


T 

S/0 

544 K 

Po 

10070 kg m~ 2 

c p 

144.87 Jkg-'K- 1 

V 

1.837 xlO" 7 mV 

D 

2.7 xlO -9 mV 1 

Co 

1 at% 

Pr 

1.25xl0“ 4 K~ l 

Pc 

-0.3049 {volume fraction )' 

K 

0.29 at%/at% 

m 

2.32 K/at% 

So 

9.8 ms' 2 

H 

5.25 xlO 4 Jkg 

k 

12.4 Wm~'K~' 


6.5 Wm~ l K~ l 

k 

amp 

2.01 

K 

8.5 xlO” 6 mV 1 

Xcry 

4.5 xlO" 6 mV 1 

amp 

8.7 xlO -7 mV 1 

S n% 

0.8 

amp 

System parameters 

W; 

6.0x10 3 m 

K 

1.0 xlO" 2 m 

L 

4.2 xlO" 2 m 


L h 4.95x1 (T 3 ot 

L e 4.95 xlO" 3 m 

L g 3.21 xlO" 3 m 

dT/dx 2.0 xlO 4 Km~ ] 

T h 913 K 
T c 323 K 
V g -3.34x10 ~ 6 ms' ] 

B 0 0.5 T 

Scale parameters 

U 0 6.12xl0“ 5 ms~ ] 

AT 650 K 
R 0 3x 10~ 3 m 
t 0 49.0 s 

Nondimensional parameters 

Pr 0.0215 
Sc 68.0 
Gr T 6.37 xlO 5 
Gr s - 1.957 xlO 4 

Ste 0.558 
Ha 38.7 
Ra 0.01859 

T m 0.34 
V p -0.055 
T' c 0.497 



Figure Captions 

Figurel. Schematic representation of the 2-D model for the Bridgman-Stockbarger Sn-doped Bi 
single crystal growth system in microgravity. 

Figure2. Dependency of interface position on concentration striation along the interface: (a) the 
interface positions, the temperature and concentration distributions along the growth interface of 
the steady cases with and without concentration effect respectively and (b) their differences. 

Figure3. The development of the flow field caused by a single frequency g-jitter in one time 

harmonic cycle in the absence of magnetic field: (a )t = 3.876, U mix = 0.6161 

(b)/ = 3.927, U max = 0.6922 (c)t = 3.978, U max = 0.6161 (d) t = 4.029, U m3X = 0.6920 . 


Figure4. The development of the concentration distribution caused by a single frequency g-jitter 
in one time harmonic cycle in the absence of magnetic field: (a) r = 3.876 (b)t = 3.927 
(c)t = 3.978 (d) t = 4.029. From left to right, C max = 3.52 reduced to C = 1.06, AC = 0.35 

Figure5. The development of (a) the interface positions, (b) the temperature and (c) 
concentration distributions along the growth interface caused by a single frequency g-jitter in one 
time harmonic cycle in the absence of magnetic field. 

Figure6. The development of the flow field caused by a single frequency g-jitter in one time 

harmonic cycle in the presence of magnetic field: (a )t = 3.876, U max = 0.1603 

(b)f = 3.927, U mxx = 0.5483 (c)f = 3.978, U m = 0.1603 (d) t = 4.029, U w = 0.5483 . 


Figure7. The development of the concentration distribution caused by a single frequency g-jitter 
in of one time harmonic cycle in the presence of magnetic field: (a) t = 3.876 (b) t = 3.927 
(c)t = 3.978 (d) t = 4.029 . From left to right, C max = 3.50 reduced to C = 1.04 , AC = 0.35 

Figures. The development of (a) interface positions, (b) the temperature and (c) concentration 
distributions along the growth interface caused by a single frequency g-jitter in of one time 
harmonic cycle in the presence of magnetic field. 

Figuic9. The effects of the magnotic field on the velocity and concentration oscillation caused by 
a single frequency g-jitter at selected points in the melt near the growth front (a) x=3.6, y=0.5 (b) 
x=3.6, y=0.5 (c) x=3.6, y=0.1 and x=3.6,y=1.3. 


Figurel 0. The effects of the magnetic field on the velocity and concentration oscillation caused 
by real g-jitter data at selected points in the melt near the growth front (a) x=3.6, y=0.5 (b) x=3.6, 
y=0.5 (c) x=3 .6, y=0.1 and x=3.6,y=1.3. 


Figurel 1. The development of (a) growth interface positions and (b) the concentration 
distributions along the growth interface caused by real g-jitter data without and with magnetic 
field at the selected time. 



Figurel2. The time variations of the concentration uniformity along the growth interface, which 
measured by (c(t) max - C(f) min )/C(t) flVSrage with and without the applied magnetic field 
respectively. 
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